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1 Introduction 

To model the sources of measured magnetic fields, the 
relationship of the magnetic sensors to the head must 
be known [2], This is usually accomplished by feed¬ 
ing currents through coils attached to the surface of 
the head [4], The locations of the coils with respect 
to the sensor array are then computed on the basis of 
the measured signals. 

When evoked responses are studied, several responses 
usually need to be averaged. Usually the position of 
the head is assumed to be fixed during the experi¬ 
ment, but if the head has moved or data from differ¬ 
ent measurement sessions are combined, the combi¬ 
nation of different head positions blurs the data spa¬ 
tially and adds a location bias. Therefore, the move¬ 
ments should be either eliminated or taken into ac¬ 
count in the analysis. 

Subjects are typically asked to keep their head still 
during the measurements. This may be difficult even 
for cooperative subjects if the measurement is long 
or subject’s task includes movements. Using indi¬ 
vidual bite-bars helps avoiding head movements [6], 
but some subjects find this rather uncomfortable and 
it can not be applied in experiments requiring verbal 
responses. 

If the effects of measured head movements can be 
corrected for during data analysis, measuring MEG 
of less cooperative subjects, such as children, will be 
easier and there will be no additional restrictions on 
subject’s tasks. In this study we compare three differ¬ 
ent approaches for this kind of head movement cor¬ 
rection. 

2 Methods 

2.1 Minimum norm estimate with singular 
value truncation 

Let us consider the minimum norm estimate [3] regu¬ 
larized with singular value truncation: 

min ||q|| subject to 

T(a)G(a:)q = T(a:)b, (1) 


where q is a ^-dimensional vector representing cur¬ 
rent strengths, G(x) is the rn, x k gain matrix with 
head position x, and b is the measured data. Matrix 
G(x) has the singular value decomposition 

G(x).= U'{ig)A(a:)V^)*, (2) 

and T(:r) consists of n first singular vectors of G(x): 

T(x) = U n (x) t . (3) 

With the assumption that the head position has been 
xo during all epochs, the minimum norm estimate is 

q = V„(a:o)A ?l (a:o) + U n (a:o)* (4) 

To correct the head position MNEs can be calculated 
separately for each epoch and averaged afterwards 
leading to the MNE-corrected estimate: 

Qmne = ^E V »^)An(s) + U n (^b. (5) 

This may be a rather time-consuming process if each 
epoch has a different head position. Fortunately, a 
faster method can be constructed by modifying the 
regularization. Since 

T(a) G(a) = A n (x) V n ($){ (6) 

the regularization can be thought of as a projection 
of the data to that generated by a n-dimensional sub¬ 
space L(x) of the source space: 

L(i)={V„(i)c:ceR"}. (7) 

The source subspace L(x) depends, in principle, on 
the head position. However, if the estimate is reg¬ 
ularized with the source subspace of a representa¬ 
tive head position xq, the small, n-dimensional source 
space L{xq) can be averaged instead of the large, k- 
dimensional source space. 

For different head positions the representative source 
space L(xo) should still explain most of the data: 

b ps G(a;)V n (a:o) c (8) 





and c can be estimated as 

c = (G(z) V„(zo)) + b. (9) 

The fields in the representative head position are 
b* = G(zo)V n (a:o) (G(a) V„(zo)) + b, (10) 
and the regularized field is correspondingly 

T(z 0 )b* = U^zo^G^V^o) 

(G(a;) V n (a; 0 )) + b 
= A n (a; 0 ) (G(a;) V„(a;o)) + b.(ll) 

The regularizing matrix transforming the fields gen¬ 
erated in head position x to those generated in head 
position xq is thus 

T(*r = An(igo) (G(ar) V n (z 0 )) + , (12) 

and the estimated source distribution using the signal- 
space correction is 

4s = V n 0r o )A„(z o ) + ^-— X^T^b- ( 13 ) 

While equation 5 requires for each head position the 
singular value decomposition of the m x k gain ma¬ 
trix, the equation requires only a pseudoinverse of a 
smaller n x n matrix, and is thus computationally 
much more efficient. 

A third correction method is to use the ordinary aver¬ 
age and simulate the head movements in the forward 
calculation. That is, both the matrix G(a:) and the 
data vector b are replaced with their averages, lead¬ 
ing to forward calculation corrected estimate 

q fw = V n A+U* — V b, (14) 

n ave ^ 

where UAV* is the singular value decomposition of 
the average gain matrix ]T G(.x). The blurring 
of fields with different head positions remains, but the 
bias caused by the head movements is removed with¬ 
out adding the bias caused by the regularization. 
When the depth weighted minimum norm estimate 

min ||W(a:) q|| subject to 
T(s)G(a:)q = T(a:)b, (15) 

is applied, where W(rr) is a diagonal weighting ma¬ 
trix compensating for different sensitivities at differ¬ 
ent source locations, the respective formulas can be 
derived using the change of variables q' = W q. 


For example, the regularized transformation (Eq. 12) 
changes to 

T(x)* = A„ (xo) (G(x)W(x)- 1 V n (x 0 )) + , (16) 

where A and V are calculated from the normalized 
gain matrix GW -1 . The respective estimated source 
distribution (Eq. 2.1) is 

4s = W(a; 0 ) _1 V n (a; 0 )A n (xq) + ]T T(z)*b. 

(17) 

2.2 Realigning magnetic fields 

The minimum norm estimates can also be used to 
transform the measured fields to fields representing 
a standard space [5]. The basic idea is to calculate the 
magnetic field that would be generated in the standard 
head position by the estimated currents. The trans¬ 
formed field using MNE-correction (Eq. 5) is 

b M NE = G(a:o)^X) V «( a: ) A «( a: ) +U «( a: ) tb - 

(18) 

The transformed field using the signal-space correc¬ 
tion (Eq. 2.1) can be calculated as 

b ss = V n (x 0 )^E T ( a; )* b ( 19 ) 

and the forward calculated correction (Eq. 14) as 

b fw = G(x 0 ) V n A+U^ ^ E b - ( 2 °) 

2.3 Multidipole modeling 

The realigned magnetic fields can then be used for 
dipole modeling. Since the signal-space correction 
yields fields lying in the subspace spanned by the 
columns of U n (io)> the corresponding projection op¬ 
erator must be applied to the dipole model fields when 
the least-squares error function is calculated. 

Some parts of the field must, however, be extrapo¬ 
lated, and the accuracy of the transformation will vary 
between different virtual magnetometers. This should 
be taken into account in the inverse problem calcula¬ 
tions. 

The forward calculated correction can be applied di¬ 
rectly to dipole modeling without transforming the 
fields simply by replacing the dipole forward calcula¬ 
tion with the average field of different head positions. 



2.4 Simulations 

The methods were evaluated by using simulated mea¬ 
surements with a neuromagnetometer with 122 pla¬ 
nar gradiometers Neuromag-122™ [1]. The simu¬ 
lated head movements consisted of a sinusoidal verti¬ 
cal movement with an amplitude of 20 mm. Since the 
order of different head positions does not contribute to 
any of the proposed methods, this closely resembles 
a possible real movement where the subject’s head 
gradually moves downward during an experiment. 
The simulated data consisted of signals produced by 
single dipoles in 400 places and two orientations, 
active one at a time. The conductivity model was 
spherically symmetrical and the dipole locations were 
50 mm from the center. A 3-dimensional cubical lat¬ 
tice with 10 mm distance between neighboring points 
was used as a source space in the minimum norm es¬ 
timate. 70 singular values were used after the regular¬ 
ization. 

The accuracies of the correction methods were stud¬ 
ied without additional noise. Assuming the measure¬ 
ment noise is white, the linear correction methods in¬ 
creases the noise level in proportion to the Frobenius 
norm of the transformation matrix, and this was used 
to evaluate the noise sensitivity of different methods. 
The head position for the non-corrected simulations, 
as well as the reference head position into which the 
different alignments were made, was the average head 
position, which is also the best ra nk -one approxima¬ 
tion for the different methods. In real measurements 
the head position is typically measured before the ex¬ 
periment, which leads to a bigger bias, if the head 
moves. Therefore the use of the average head posi¬ 
tion can be viewed as an approximate head position 
correction, which is not, however, as accurate as the 
proposed correction methods. 

3 Results 

The accuracy of the minimum norm estimates using 
different corrections were evaluated using the loca¬ 
tions of the maximal signals. Even though the MNE 
is an unbiased estimator of the current strength, the lo¬ 
cation of the maximal value of the regularized MNE is 
not a perfect estimate of the location of a focal source, 
and even with a correct head position the root-mean- 
square (RMS) location error was 17 mm. To over¬ 
come this source of errors the maximum location of 
the corrected value from movement simulations was 
compared with that of the MNE in the reference head 


position (Fig. 1). 

The RMS difference over different simulated dipoles 
was over 16 mm when no correction was applied 
and dropped below half of that using any correction 
method. The noise level (Frobenius norm of the trans¬ 
formation matrix) was close to original value using 
the forward correction, and largest using the signal- 
space correction. 
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Figure 1: The accuracy of the MNE maxima and the 
noise level with different correction methods 

The accuracy of the magnetic field alignment was 
evaluated similarly. The accuracy of the correction 
was measured as the relative RMS difference between 
the realigned magnetic field and the simulated field 
with the reference head position. Due to the regu¬ 
larization, the realignment methods only estimate the 
fields of the signal subspace. Therefore the noise level 
of each method was measured as the Frobenius norm 
of the transformation matrix related to the signal sub¬ 
space in the reference head position. 

The RMS difference of the fields was 19 % when no 
correction was applied, 11 % using the signal-space 
correction, and about 2 % when the MNE or forward 
correction was used (Fig. 2). The noise levels were 
similar to the values in the MNEs. 
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Figure 2: The accuracy of the realigned magnetic 
fields and the noise level in the signal space 

The accuracy of dipole modeling was evaluated by 







least-squares-fitting one dipole to each simulated (and 
possibly corrected) field. Even the uncorrected dipole 
locations were rather close to the correct values, with 
an RMS error slightly above 5 mm. The MNE and 
signal-space corrections decreased the error below 
1 mm, and the forward correction, having the correct 
system matrix, gave exactly the right locations (Fig. 
3). 



Figure 3: The accuracy of the dipole fitting in the 
noiseless case and the noise level of corresponding 
magnetic fields 


The dipole noise level was measured as the noise level 
of the corresponding realigned magnetic field. For 
forward correction no field alignment is needed, an 
therefore the amplitude noise level not increased. The 
effect of the different correction methods to the noise 
sensitivity of the dipole location parameters modeling 
was not studied, but it is likely that the small noise 
level of the forward correction is balanced to some 
extent with somewhat bigger sensitivity of the dipole 
location parameters. 

The calculation time for the MNE correction was 
largest, since that included handling all the evoked re¬ 
sponses separately and a separate pseudo-inverses for 
system matrices in each head position. The signal- 
space correction still needs to go through individual 
epochs, but only the pseudo-inverse of the system 
matrix in the reference position is needed, and was 
thus clearly faster. The forward correction was cal¬ 
culated from the normal average and did not increase 
the computation time at this stage. The forward cal¬ 
culations are slower, but the total time is negligible 
compared to that of the other correction methods. 


timates the errors caused by head movements were 
smaller than the contribution of other error sources in 
typical MEG studies. The forward correction did not 
significantly increase the time of the calculations or 
noise level of the estimates, and is thus a suitable way 
of correcting for the head movements in MEG exper¬ 
iments. 


Acknowledgements 

Supported by the Ministry of Education. 


References 

[1] A. I. Ahonen, M. S. Hamalainen, M. J. Kajola, 
J. E. T. Knuutila, R R Laine, O. V. Lounasmaa, 
L. T. Parkkonen, J. T. Simola, and C. D. Tesche. 
A 122-channel SQUID instrument for investigat¬ 
ing the magnetic signals from the human brain. 
Physica Scripta, T49:198-205, 1993. 

[2] M. Hamalainen, R. Hari, R. J. Ilmoniemi, 
J. Knuutila, and O. V. Lounasmaa. Magnetoen¬ 
cephalography — theory, instrumentation, and 
applications to noninvasive studies of the work¬ 
ing human brain. Rev. Mod. Phys., 65(2) :413- 
497, 1993. 

[3] M. S. Hamalainen and R. J. Ilmoniemi. Interpret¬ 
ing magnetic fields of the brain: minimum norm 
estimates. Med. & Biol. Eng. & Comput., 32:35- 
42, 1994. 

[4] J. Knuutila, A. I. Ahonen, M. S. Hamalainen, 
R. J. Ilmoniemi, and M. J. Kajola. Design con¬ 
siderations for multichannel SQUID magnetome¬ 
ters. In H. D. Hahlbohm and H. Liibbig, editors, 
SQUID’85: Superconducting Quantum Interfer¬ 
ence Devices and their Applications, pages 939- 
944, Berlin, 1985. Walter de Gruyter. 

[5] J. Numminen, S. Ahlfors, R. Ilmoniemi, J. Mon- 
tonen, and J. Nenonen. Transformation of mul¬ 
tichannel magnetocardiographic signals to stan¬ 
dard grid form. IEEE Trans. Biomed. Eng., 
42(l):72-77, 1995. 


4 Discussion 

The results show that using any of the proposed meth¬ 
ods the errors in dipolar or distributed current esti¬ 
mates can be decreased considerably. In corrected es¬ 


[6] K. D. Singh, I. E. Holliday, P. L. Furlong, 
and G. F. A. Harding. Evaluation of MRI- 
MEG/EEG co-registration strategies using Monte 
Carlo simulation. Electroenceph. Clin. Neuro¬ 
physiol., 102:81-85, 1997. 




